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We derive a general master equation relating the gravitational- wave observables r and Q.^^ (f); or 
the observables i}^" (fi) and r2Q™(/2). Here r is the so-called "tensor-to-scalar ratio," which is con- 
strained by cosmic-microwave-background (CMB) experiments; and ^l^" (f) is the energy spectrum 
of primordial gravitational-waves, which is constrained e.g. by pulsar-timing (PT) measurements, 
laser-interferometer (LI) experiments, and the standard Big Bang Nucleosynthesis (BBN) bound. 
Differentiating the master equation yields a new expression for the tilt d In n^™ (/)/d In / of the 
present-day gravitational- wave spectrum. The relationship between r and (f) depends sensi- 
tively on the uncertain physics of the early universe, and we show that this uncertainty may be 
encapsulated (in a model-independent way) by two quantities: w{f) and ht{f), where rit(/) is a 
certain logarithmic average over n^ik) (the primordial tensor spectral index); and w(/) is a certain 
logarithmic average over i&(a) (the effective equation-of-state parameter in the early universe, af- 
ter horizon re-entry). Here the effective equation-of-state parameter w{a) is a combination of the 
ordinary equation-of-state parameter w{a) and the bulk viscosity C,{a). Thus, by comparing observa- 
tional constraints on r and ^^{f), one obtains (remarkably tight) constraints in the {«)(/), ^^t(/)} 
plane. In particular, this is the best way to constrain (or detect) the presence of a "stiff" energy 
component (with w > 1/3) in the early universe, prior to BBN. (The discovery of such a compo- 
nent would be no more surprising than the discovery of a tiny cosmological constant at late times!) 
Finally, although most of our analysis does not assume inflation, we point out that if CMB experi- 
ments detect a non-zero value for r, then we will immediately obtain (as a free by-product) a new 
upper bound w < 0.55 on the logarithmically averaged effective equation-of-state parameter during 
the "primordial dark age" between the end of inflation and the start of BBN. 

PACS numbers: 



I. INTRODUCTION 

A variety of different experiments (some already oper- 
ating, others in various stages of development) are hop- 
ing to detect gravitational waves (tensor perturbations) 
from the early universe. In particular, at long wave- 
lengths, cosmic- microwave-background (CMB) experi- 
ments [il, 0, i, i, i, i, 0, 0, 0, [H [nl, [H, [li will measure 
(or tightly constrain) the so-called tensor-to-scalar ratio 
r by searching for its characteristic "B-mode" imprint 
in the CMB polarization anisotropy [3, [H, [l^- And 
on shorter wavelengths, various techniques — including 
pulsar-timing (PT) ITI. [TsI. [l9l and laser-interferometer 
(LI) experiments [MlEriilia, [ 



will 



sure or constrain the present-day gravitational-wave en- 
ergy spectrum r2§™(/). 

The coming decade is likely to see exciting progress 
in this area. At the lowest frequencies, CMB polariza- 
tion experiments will either dete ct g ravitational waves 
from inflation [13, [H, [H, [10, [Ml, [sj [H, [13, [H, [H] , or 
else rule out the simplest (and arguably the most com- 
pelling) inflationary models |36|. At intermediate fre- 
quencies, pulsar timing arrat;s[18l [l9| will reach far be- 
yond the gravitational-wave sensitivity of individual pul- 
sars. And at high frequencies, the sensitivity of ground- 



based gravitational-wave detectors (and also the space- 
based mission LISA, if it is launched) will surpass the 
so-called "standard Big Bang Nucleosynthesis (sBBN) 
bound" by several orders of magnitude, and thus place 
genuinely new constraints on the primordial gravitational 
wave signal at high frequencies. 

Since primordial gravitational waves provide a rare and 
precious window onto the extremely-high-energy physics 
of the infant universe, it is essential to think carefully 
about the information that they carry. 

In this paper we will present a general (yet rather 
simple) master equation ([T|) connecting the long- 
wavelength observable r to the short-wavelength observ- 
able ri§^(/). The goal is to clarify the relationship be- 
tween gravitational-wave constraints at different wave- 
lengths, and to highlight the important and unique in- 
formation about the early universe that is encoded in 
this relationship. 

What exactly do we learn, in general, by comparing 
long-wavelength constraints on r and shorter-wavelength 
constraints on {f)l From the master equation ([1]), 
we will see that this type of comparison should be in- 
terpreted as primarily constraining two quantities, w(/) 
and fitif), which encode information about the early uni- 
verse in a model-independent way. These two quantities 
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are defined by Eqs. ([7]) ~ ([9]); and explained in detail 
in Sec. IIIII For now, let us briefly discuss their phys- 
ical meaning: is the logarithmic average (over a 
certain range of comoving wavenumber k) of the primor- 
dial tensor spectral index n^{k); and w{f) is the loga- 
rithmic average (over a certain range of the cosmological 
scale factor a) of the effective equation-of-state param- 
eter w{a) in the early universe (after horizon re-entry). 
Here the effective equation-of-state parameter w{a) is a 
combination of the ordinary equation-of-state parameter 
w{a) and the bulk viscosity ((a): see Eq. 

A key advantage of our current formulation in general 
(and of the variables w{f) and fitif), in particular) is 
that w{a), C(a) and n((fc) may be arbitrary functions of 
a and A:, respectively. So, in particular, we will not take 
w or Uf. to be constant (or piecewise constant), as is often 
assumed in analytical treatments of primordial gravita- 
tional waves. The point is that, when deriving Eq. p]). 
the quantities w{f) and fitif) naturally arise as the most 
direct and general encapsulation of the uncertain early- 
universe physics that enters into the relationship between 
r and 17r(/). 

As an application, we will stress that comparison r and 
^o^(/) provides the most powerful way to constrain the 
equation-of-state parameter w{a) during the "primordial 
dark age." Here we use the phrase "primordial dark age" 
to refer to the epoch separating the end of inflation from 
the start of Big Bang Nucleosynthesis (BBN). Note that, 
on a logarithmic scale, this primordial dark age spans a 
large fraction of cosmic history: the energy scale of BBN 
is ^ 10~^ GeV, while the energy scale at the end of infla- 
tion may exceed 10^^ GeV. And yet, although there is a 
standard theoretical picture of how the universe behaves 
during this early epoch, we currently have essentially no 
direct observational constraints. 

In fact, there are several reasons to be nervous about 
one of the key (implicit) assumptions in the standard 
picture of the primordial dark age: namely, the assump- 
tion that the equation-of-state satisfies w < 1/3. The 
first reason to worry is rather general: since the en- 
ergy density of a cosmological matter component scales 
as p oc a~^^^^'^\ components with lower w dilute more 
slowly. Thus, just as an exotic component with w suffi- 
ciently low will tend to dominate the cosmic energy bud- 
get at sufficiently late times (think of "dark energy" with 
w < —1/3), an exotic component with w sufficiently high 
(call it "stiff energy" with w > -1-1/3) will tend to domi- 
nate the cosmic energy budget at sufficiently early times 
(see Fig.[T]). Indeed, as we look backward past BBN, the 
primordial dark age provides a huge window in which a 
stiff energy component might overtake radiation as the 
dominant component in the cosmic energy budget, with- 
out coming into conflict with any current observational 
constraints. It is also worth noting that there are per- 
fectly sensible energy components with w > 1/3 which 
might be present in the early universe. For example, a 
homogeneous scalar field (/<(t) with vanishing (or negligi- 
ble) potential energy V{(j)) = has w = 1; and, in fact. 
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FIG. 1: How the components of the cosmological energy bud- 
get scale with cosmological expansion: "stiff energy" (solid 
purple line), radiation (long-dashed red line), matter (dotted 
blue line), and dark energy (dot-dashed green line). Compo- 
nents with higher w tend to dominate at earlier times. Our 
universe may be dominated by a "stiff energy" component 
(with w > 1/3) prior to Big Bang Nucleosynthesis (but after 
inflation) . 

supergravity and string theory seem to naturally pre- 
dict many (embarrassingly many!) scalar moduli fields 
with precisely this property. Furthermore, various au- 
thors [33, Ha, [3^, lio, |4l|, have considered inflation 
models in which the inflaton field itself experiences a 
period of free {w = 1) evolution at the end of infla- 
tion; or some other equation-of-state stiffer than radi- 
ation [n, H, m, H gzi, m . 

"Stiff {w > 1/3) energy" in the early universe may 
seem like an exotic possibility. But would the discovery of 
"stiff energy" at early times be any more surprising than 
our apparent discovery of "dark energy" at late times? 
One lesson that we have learned from dark matter and 
dark energy is that the universe has an unmistakable pen- 
chant for new and unexpected energy components; and 
it is important to check for these components observa- 
tionally, if possible, rather than simply assuming that 
they are not there. We will stress that the comparison of 
constraints on r and CIq"" (f) provides the best means for 
carrying out such a check. 

One of the most important results in this paper comes 
from considering the relationship between the CMB con- 
straint on r and the sBBN constraint on ilf"{f). If CMB 
polarization experiments succeed in detecting a non-zero 
value for the primordial tensor-to-scalar ratio r, this will 
be widely interpreted as providing evidence for inflation. 
But we show that, if these primordial tensor fluctua- 
tions are really generated by inflation, then (in combi- 
nation with the current sBBN constraint on nf"{f)), 
this will also imply an immediate and important sup- 
plementary result: namely a remarkably tight bound in 
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the plane. This bound in the {w{f),'h^{f)} 

plane is derived in Sec. IVIIi and shown in Fig. [5l 

If CMB polarization experiments detect a non-zero 
value for r, then the bound depicted in Fig. [5] will be 
a qualitatively new piece of model-independent informa- 
tion about the early universe — which is very exciting, 
since such information is notoriously hard to obtain! One 
way to look at the bound is as follows: If we assume that 
the bulk viscosity ((a) is negligible after inflation, and 
also that the primordial tensor power spectrum A^(fc) is 
nearly flat (which is a prediction of inflation), then we 
obtain an upper bound {w) < 0.55 on the logarithmic 
average of the equation-of-state parameter w{a) during 
the primordial dark age separating the end of inflation 
from the BBN epoch. 

It is important to clarify the range of validity of our 
analysis. When we use f2g^(/) in this paper, we are refer- 
ring only to primordial gravitational waves — and, more 
specifically, only to those gravitational waves that were 
generated well before the corresponding comoving wave- 
length "entered the Hubble horizon" (i.e. became shorter 
than the instantaneous Hubble length). Apart from 
this restriction, the results are quite general, and make 
no assumptions about the physical mechanism respon- 
sible for generating the gravitational waves. For exam- 
ple, our analysis applies to the primordial gravitational- 
wave spectrum generated during inflation; and it ap- 
plies equally well to the primordial gravitational-wave 
spectra generated by the "pre-big-bang" [iol . [soj and 
"cyclic/ekpyrotic" [Hll, [H, HI] alternatives to inflationary 
cosmology; and, although all of the previous three exam- 
ples (inflationary, pre-big-bang, and ekpyrotic/cyclic cos- 
mology) generate primordial gravitational waves through 
the cosmological amplification of quantum fluctuations, 
our analysis would also apply to models that generate a 
primordial gravitational-waves via some completely dif- 
ferent mechanism (as long as they are generated prior to 
horizon entry). ^ On the other hand, our analysis does not 
apply, e.g. to the gravitational- wave spectrum produced 
by the incoherent superposition of signals from merging 
binary stars (53 | , or by a hypothetical period of preheat- 
ing after inflation [5a|, or by bubble collisions after a 
cosmological phase transition [5^ — since all of these 
production mechanisms result in gravitational waves that 
are shorter than the instantaneous Hubble length at the 
time when they are generated. 

This paper is organized as follows. Sec. HIl introduces 
some notation. In Sec. IIIII we present and explain the 
master equation ^ which relates r and ^lQ^{f), and will 
serve as the basis for most of our analysis. In Sec. llVlwe 



^ A caveat is that the derivation of Eq. 1131 1 specifically applies 
to standard inflation 128^ (and not to pre-big-bang 1491 or ekpy- 
rotic/cyclic ISll . |52|| models). But this is a very mild caveat, 
since Eq. JTSj is used only in Sec. IVIII which deals with models 
that produce a detectable value for r (which pre-big-bang and 
ekpyrotic/cyclic models do not ISOl l5ll Issll ). 



use the master equation ([Ij to derive two simple results. 
The first result is Eq. (fTS)) which expresses the relation- 
ship between f2Q^(/i) and ^^^^(/j) — that is, two differ- 
ent short-wavelength constraints {e.g. from LIGO and 
LISA) at two different frequencies fi and /2 . The second 
result is Eq. p9p . which significantly generalizes previous 
expressions for the tilt dlnilf" (f) / d\n k of the present- 
day energy spectrum of primordial gravitational- waves. 
In Sec.|V]we analyze the implications of combining CMB 
constraints on r with LI and PT constraints on nf"{f). 
The section breaks into 4 parts, depending on whether we 
suppose that CMB and LI/FT experiments have detected 
the gravitational-wave background, or merely bounded 
it from above. The results are summarized in Figs. 2, 3 
and 4. In Sec. lVIi we discuss the following point: If LI or 
PT experiments detect an unexpectedly strong stochastic 
gravitational- wave signal flf"{f), then there is a rough 
observational consistency check that this signal should 
satisfy, if it is truly of primordial origin. In Sec. IVIII we 
analyze the constraint in the {w{f),nf{f)} plane that 
follows from combining a CMB detection of r with the 
sBBN bound on r2g™(/). As mentioned above, this con- 
straint is rather strong; and it is also quite insensitive to 
the detected value of r: see Fig. 5. Finally, we conclude 
in Sec. IVIIII Some of the key equations in the text are 
derived in appendices: in particular, Eq. ([1]) is derived in 
Appendix [XI and Eqs. and are derived in Ap- 
pendix[B] Appendix [Cl lists a few numbers that are useful 
for converting our various algebraic expressions into nu- 
merical results and plots. 



II. NOTATION 

Throughout this paper, we will often use subscripts 
to indicate the time at which a quantity is to be evalu- 
ated. For example, a quantity with subscript "0" is eval- 
uated at the present time; a quantity with subscript "eq" 
is evaluated at the moment of matter-radiation equality 
(PocT* = Poq"^); a quantity with subscript "c" is evaluated 
at the redshift Zc (defined in Sec. lIIip : and a quantity with 
subscript "fc" is evaluated when the comoving wavcnum- 
bcr k "re-enters the Hubble horizon" {i.e. crosses from 
k < aH to k > aH). 

We will also use units in which the speed of light is 
unity, c = 1. 



III. THE MASTER EQUATION 

Primordial gravitational wave measurements probe 
two basic quantities. On long wavelengths, CMB polar- 
ization experiments constrain the tensor-to-scalar ratio 
r. And on shorter wavelengths, various techniques con- 
strain the present-day gravitational-wave energy spec- 
trum Vlf^{f). In Appendix 1X1 we derive a master equa- 
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tion relating r and rig^(/). The result is: 



(1) 



As we shall see in a moment, the factor Ai is roughly 
independent of the gravitational- wave frequency /, while 
the two factors A2 and A3 are both proportional to /, so 
that nl^if) is roughly proportional to f^if^+'hif). 

Now let us carefully explain the meaning of each quan- 
tity appearing in Eq. ([T]) — namely, the gravitational- 
wave observables {r2Q*(/), r}, the factors {^1,^2, A3}, 
and the exponents {a(f),n^{f)}. 

The present-day gravitational-wave energy spectrum 



1 dpi- 
pf'dlnf 



(2) 



represents the present-day gravitational-wave energy 
density {pf") per logarithmic frequency interval, in 
units of the present-day "critical density" yo"'* = 
3H^/{8TrGj^), where i?g is the present-day value of the 
Hubble expansion rate, and is Newton's gravitational 
constant. 

The tensor-to-scalar ratio 



(^cmb) 
^TJ.(^cmb) 



(3) 



is the ratio of the primordial tensor power spectrum 
^h(^cmb) (defined in Appendix A) to the primordial 
scalar power spectrum A|j(fc^j^,^) at the CMB wavenum- 



ber k^^^. Our definition of the tensor-to-scalar ratio 
matches the convention used, e.g., by the WMAP ex- 
periment [13, \EM CAMB numerical code [soj : 
but beware that there are several alternative defini- 
tions/conventions floating around in the literature. The 
CMB wavenumber k^^^ is the comoving wavenumber at 
which CMB experiments report their constraints on A^, 
Al, and r: e.g. the WMAP experiment [13, [HI] uses 
= 0.002 Mpc-\ where ao is the present-day 
value of the cosmological scale factor. 

Next consider the 3 factors {^1,^2,^3} appearing in 
Eq. 11]). They are given by 



Ai 



_C2(fc)C3(fc)A^(fc,„,j7 
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A2 = 
A3 = 



1 



Hq J (l + z,)7i/2' 
2^/^ Ho 



Ho J (fccmb/«o) ' 



(4a) 
(4b) 

(4c) 



where 



^ Note that, in Eqs. JTJ and Il4al l. we could trade the more- 

I : 

less commonly used) observable A^(fc^j^j^) 



commonly-used observables r and A|j(fc(,mb) ^^'^ single (but 



Here Q^''^ 



P6 



™'^*/p"'* is the ratio of the present-day 



non-relativistic matter density p™^* to the present-day 
critical density /c"'*. The comoving wavenumber k is 
related to the physical frequency / through the rela- 
tion k/ao = 27r/. The cosmological scale factor a is 
related to the cosmological redshift z through the rela- 
tion ag/a = 1 + z. In particular, Zeq denotes the red- 
shift of matter-radiation equality (p™*** = Pe^^), while 
Zc denotes the highest redshift at which we know that 
the universe was radiation dominated {i.e. the redshift 
at the end of the "primordial dark age" discussed in the 
introduction). Given our present observational knowl- 
edge of the early universe, it is natural to choose Zc 
to be the redshift of BBN, z^^^^; but in the future, as 
our knowledge of the early universe improves, a differ- 
ent choice (i.e. a higher redshift z^) may become more 
appropriate. The factors g*{z) and g^siz), which mea- 
sure the effective number of relativistic degrees of free- 
dom in the universe at redshift z, are conveniently de- 
fined as follows: If p{z), s{z), and T{z) denote, respec- 
tively, the energy density, entropy density, and temper- 
ature at redshift z, then p{z) = {tt^ / 30) g^,{z)T'^{z) and 
s{z) = {2Tr^/ 45) g^g{z)T^{z). For a detailed discussion 
of the "correction factors" C2{k) and C3{k), including 
definitions and explicit expressions, see Appendix El and 
Rcf. [g^]. For now, it is enough to note that C2(fc) and 
C3{k) are both 0(1), which means that they will not 
play a very significant role in this paper (although in 
other contexts they can be interesting and important, 
see Ref. ^). 

Finally consider the two exponents, a{f) and fifif), 
that appear in Eq. ((T|). The first exponent, a{f), is given 

by 



3w{f) + l 
where u)(/) is the logarithmic average 

1 



Hf)^ 



ln(ac/afc) Ja 



_ , .da 
w[a) — 



(6) 



(7) 



of the effective equation-of-state parameter w{a) from Ofc 
(the scale factor when k = 27rao/ re-entered the horizon) 
to flc (the scale factor at redshift Zc). Here the effective 
equation-of-state parameter w(a) is given by 



■w{a) = w{a) 



87rG'jvC(a) 
Hia) ■ 



(8) 



where w{a) = p(a)/p{a) is the ordinary equation-of- 
state parameter [i.e. the ratio of the total cosmological 
pressure p{a) to the total cosmological energy density 
(5) p{a) = p'^"'(a)], H{a) is the Hubble expansion rate, and 
C(a) is the bulk viscosity of the cosmological fluid (see 
Sees. 2.11 and 15.11 in Ref. [6l|). The second exponent, 
fij(/), is given by the logarithmic average 



1 



ln(fc/fcc.nb) 



nt(fc') 



dk' 



(9) 
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of the primordial tensor tilt nf{k') over the wavenumbcr 
range fccmb < k' < k. Here the primordial tensor tilt 
nt(fc) is defined as the logarithmic slope of the primordial 
tensor power spectrum Af^{k) at comoving wavenumber 
k: 



nt{k) 



^ dlnA|(fc)^ 
c?ln fc 



(10) 



We again stress that cquation-of-statc parameter w{a) 
may have arbitrary a-dcpcndcncc, and the primordial 
tensor tilt n^^k) may have arbitrary fc-dependence. We 
do not assume that w or is constant. 

Let us clarify the sense in which rig^(/) is a "short- 
wavelength" gravitational-wave observable. We mean 
that, in the master equation ([T]), the quantity f2g^(/) rep- 
resents the present-day gravitational-wave energy spec- 
trum on scales that re-entered the Hubble horizon dur- 
ing the primordial dark age: that is, after the end of 
inflation (or whatever process produced the primordial 
gravitational- wave signal), but before the redshift Zc- In 
other words, the frequency / that appears in equation 
([T|) lies in the range: 



fc< f < fe: 



(11) 



Here fc is the present-day frequency of the comoving 
wavenumber kc = 27rao/c that re-entered the Hubble 
horizon (kc = OcHc) at redshift Zc, and f^^^ is the high- 
frequency cutoff of riQ^(/). As shown in Appendix B, fc 
is given by 



(12) 



and, if the primordial tensor spectrum is generated by 
inflation, then f^^^ is given by 



f =^ 

J end o 

ZTT 



cmb 



(13) 



where, in this equation, we have used the ab- 
breviated notation {a,(3,h^\ for the quantities 
{"(/o„d)'/3(/o„d)'"t(/ond)}. and defined 



/3(/c: 



"(/end) - "-t(/ond)- 



(14) 



For concreteness, let us give some rough numbers: if we 
take Zc = Zj^j,,-, (i.e. the redshift at which the temperature 
was T « 1 MeV), then fc = /bb„ ~ 1-8 x 10"" Hz; and if 
the primordial tensor spectrum is generated by inflation 
(with fit ~ 0), followed by a "standard" primordial dark 
age (with w « 1/3), then /^^^ « 4.5 x Iffr^l^ Hz. 

Let us emphasize once again that the derivation of 
Eq. ([T^ is the only place in this paper where we as- 
sume that the primordial gravitational wave spectrum 
was generated by inflation. Since most of the results in 
this paper do not rely on Eq. their validity does 

not rely on the correctness of inflation. Indeed, we will 



only need Eq. (fT3l) in Scc. lVIIl when we want to combine 
CMB and BBN constraints. 

It is useful to interpret the master equation ^ as fol- 
lows. From Eq. ([T]), we see that the relationship be- 
tween r and [f) is much more sensitive to the two 
quantities {'w{f),nt{f)} than it is to the three quanti- 
ties {Ai,yl2,j43} — because t«(/) and n^if) appear in 
the exponents of the huge dimensionlcss numbers A2 and 
A3. This means that, even though the numerical values 
of {^11,^2,^3} are somewhat uncertain (since, e.g., Hq 
and Q^^^ are measured with non- negligible error bars, 
and C2 and C3 are only known to be roughly equal 
to unity), these uncertainties do not significantly affect 
the constraints on w{f) and n((/) coming from Eq. ([1]), 
as we shall see in more detail below. In other words, 
we may think of {^1,^2:^3} as "known" quantities; so 
that when we measure or observationally constrain r and 
^o^(/)i Ihe master equation ^ allows us to directly in- 
fer constraints on the "unknown" quantities w(/) and 
ntif)- 



IV. TWO SIMPLE CONSEQUENCES 

Before moving on, we note two simple results that fol- 
low directly from the master equation ([1]). 

The first result is obtained by evaluating Eq. ([T]) at 
two different frequencies, /i and /2, and taking the ratio 
to get: 



»r(/i) _ g2(fci)C3(fci) 

nrif2) C2(fc2)C3(fc2) 

Here nt(/i, /2) is given by 



a(/l,/2)+n,(/i,/2) 



^t(/l,./2) = 



1 



dk 

nt{k) — , 



Hki/k2) 

where ki^2 — 27rag/i.2. And ci(/i,/2) is given by 

,3u;(/i,/2)-l 



with 



«(/l,/2) 



1 



ln(ai/a2) 



/ w{a) — , 



(15) 



(16) 



(17) 



(18) 



where and 02 are, respectively, the values of the scale 
factor when ki and ^2 re-entered the Hubble horizon. 

Thus, whereas Eq. ([1]) shows how long- wavelength 
(CMB) gravitational-wave constraints relate to shorter- 
wavelength (pulsar, laser-interferometer, and nucleosyn- 
thesis) constraints; Eq. (fT5|) explains how two shorter- 
wavelength constraints {e.g. from LIGO and LISA) relate 
to one another. 

The second result is obtained by differentiating Eq. ([T|) , 
which yields a new expression for the logarithmic tilt of 
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the present-day energy spectrum: 



din/ 



f3w{k)-l 
\3w{k) + 1 



d\nC2{k) dlnCsik) 



d In k 



din k 



(19) 



where w{k) is the value of the effective equation-of-state 
parameter [see Eq. (jS])] when the comoving wavenumber 
k ~ Sttoq/ re-enters the Hubble horizon (fc = aH). 

Note that Eq. (|19p generalizes earlier expressions 
0, 113, m for d\nnl^{f)/d\nf, in the sense that it 
includes the corrections arising from the following 3 phys- 
ical effects, if they are present at the moment when 
the comoving wavenumber k is re-entering the Hubble 
horizon (fc = aH) in the early universe: (i) first, the 
term involving w{k) incorporates the correction due non- 
neghgiblc bulk viscosity ( (see Eq. (|5])); (ii) second, the 
term involving C2(fc) is the correction arising from time- 
variation of the effective equation-of-state parameter w; 
and (iii) the term involving C3(fc) is the correction due 
to non- negligible tensor anisotropic stress tt^ . Again, see 
Appendix A and Ref. (g^I for more details on the correc- 
tion factors C2{k) and Cs{k). 

Furthermore, if the primordial gravitational-wave 
spectrum is produced by the amplification of vac- 
uum fluctuations as the mode k "exits the Hub- 
ble horizon" in the early universe (as in inflationary, 
cyclic/ekpyrotic, and pre-big-bang cosmological mod- 
els), and the equation-of-state parameter is varying suf- 
ficiently slowly as k exits the horizon, then n^{k) is given 

by 



n^{k) 



l-«^cxit(fc) 



1 + 3u;,,h(^) 



(20) 



(see Eq. (38) in Ref. [62]), where 'w^^i^{k) is the equation- 
of-state parameter, evaluated at the moment when k exits 
the Hubble horizon. Note that Eq. (|20[) applies equally 
well (i) to expanding models (like inflation, where the 
modes exit the Hubble horizon while the universe is ex- 
panding with w < —1/3); and (ii) to contracting models 
(like the pre-big-bang or cyclic/ekpyrotic models, where 
the modes exit the horizon while the universe is contract- 
ing with w > —1/3). 



V. CMB + LI/PT CONSTRAINTS 

In this section, we explore some of the implications 
of the master equation ([T]), focusing on the relation- 
ship between cosmic-microwave-background (CMB) po- 
larization experiments at long wavelengths and laser- 
interferometer (LI) and pulsar-timing (PT) experiments 
at shorter wavelengths. The discussion naturally breaks 
into 2x2 = 4 cases, depending on: (i) whether or not 
CMB polarization experiments have successfully detected 
r, and (ii) whether or not LI or PT experiments have 
successfully detected 17g^(/). We number these cases as 
shown in Table [H and consider each case in turn. 





LI/PT non-detection 


LI/PT detection 


CMB non-detection 


Case 1 


Case 3 


CMB detection 


Case 4 


Case 2 



TABLE I: The analysis in Sec. |V] breaks into 4 cases, de- 
pending on: (i) whether or not cosmic-microwave-background 
(CMB) experiments have already detected a non-zero value 
for r; and (ii) whether or not laser-interferometer (LI) or 
pulsar-timing (PT) experiments have already detected a non- 
(./)• 



zero value for fig" 



A. Case 1: neither r nor Q,^^{f) is detected 

First suppose that CMB experiments have not yet de- 
tected r; and that LI/PT experiments have not yet de- 
tected r2Q™(/). Of course, this is the current situation in 
2007. 

CMB observations provide an upper bound r < r^^^. 
Currently r^^^ « 0.5 [sHi- It is often claimed that 
this long-wavelength bound implies an upper bound on 
^o"{f) a-t shorter wavelengths. Let us examine this 
claim. 

In fact, from the master equation ([T]), we see that the 
upper bound is: 



where 



A, — 2 ( '^^max 1 
max I o " II 



(21) 



(22) 



In other words, in order to infer an upper bound on 
^o^(/) from the CMB upper bound r < r^^^, we need 
to assume two additional bounds: w{f) < Wj^ax ^^'^ 
fi'tif) ^ "-tmax- these two additional bounds are 

theoretical speculations about the early universe — not 
observational facts — so they should make us nervous. 
Furthermore, since A2 and A3 are huge dimensionless 
numbers, we see that the upper bound on Oq^(/) is very 
sensitive to the assumed values for w^^^ and fif ^^x- 

Now let use consider the most reasonable assumptions 
about and given our current theoretical un- 

derstanding of the early universe. 

What is the most reasonable assumption for w^^^^l 
First note that, if we assume the bulk viscosity C(a) is 
non-negative (as required by the second law of thermody- 
namics), and we assume that the equation-of-state w{a) 
satisfies the upper bound w{a) < Wmaxi then Eqs. ([7|) 
and ([8]) imply that w(/) satisfies the same upper bound: 
""^max = ''"max- Next uote that a fluid of massless (or 
extremely-relativistic) non-interacting particles satisfies 
If = 1/3; and if we give some of these particles finite 
masses, or finite interactions, this tends to decrease w 
below 1/3 (see Refs. [gO, [S^)- And in standard reheat- 
ing/preheating after inflation, one also typically finds 
w < 1/3 [131 ■ For these reasons, and others, 



1/3 



(23) 
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FIG. 2: This figure relates to "Case 1," discussed in Sec. IV Al 
The curves show the upper bound on Q,f^{f), over the range 
/cmb < / < /end) various assumcd values of i&niax and 
^t.max- The 4 solid black curves correspond (from bottom to 
top) to Wmax = {1/3, 0.4, 0.5, 0.6} and hf^^^^ — 0. The 4 dot- 
ted red curves show the same thing, but now with n^^niax ~ 
— ''max/8. The 4 dashed blue curves correspond (from bot- 
tom to top) to nt,inax = {0.1,0.2,0.3,0.4} and il^max = 1/3. 
Note that frequencies below = fy^^^^ ~ 10^^^ Hz re-entered 
the Hubble horizon after BBN, and hence are unaffected by 
assumptions about w during the primordial dark age. The 
current and future experimental constraints shown in the fig- 
ure are discussed in the text, at the end of Sec. IV Al 



is probably the best guess. But, as argued in Sec. [H 
there are also perfectly reasonable matter components 
with w > 1/3, and there are even reasons to suspect 
that these components might gcncrically be important 
at sufficiently early times (see Fig.[T]). Given our current 
understanding of the early universe, w^jax = 1/3 is a good 
guess — but it is only a guess, and should be checked 
exp er iment ally. 

What is the most reasonable assumption for ^j^x? 
First note that, if we assume that the primordial ten- 
sor tilt n^{k) satisfies the upper bound n^{k) < j^^^^, 
then Eq. ^ implies that rij(/) satisfies the same up- 



per bound: hf. 



If we assume that the pri- 



mordial gravitational-wave spectrum was generated by 
inflation, then the primordial tensor tilt is given by the 
well-known formula nj(fc) = — 2e(fc), where e{k) refers to 
the value of the parameter e{k) = (3/2)(l -I- Woxit(fc)) = 
— (i(lniJ)/(i(lna)|^ when the mode k leaves the Hub- 
ble horizon {k = aH) during inflation. Then, as long 
as the stress-energy tensor T^^ during inflation satis- 



fies the so-called "weak energy condition" (which, as its 
name suggests, is a very mild assumption, corresponding 
to w > —1 in a Friedmann-Lemaitre- Robertson- Walker 
(FLRW) universe), we can infer nj(fc) < 0. For these 
reasons and others. 







(24) 



is probably the best guess. Note that this conclusion is 
rather general within the context of inflation, in the sense 
that we have not made reference to scalar fields, or any 
other details of the (currently unknown) matter content 
driving inflation. Indeed, the conclusion should be valid 
as long as the following two conditions hold: (i) gravity 
may be described (at least effectively) by 4-dimensional 
general relativity during inflation; and (ii) e is ^ 1 and 
slowly varying during inflation. Both of these conditions 
are indeed satisfied by most viable inflationary models 
that have been considered (single-field, multi-field, . . . ), 
although there arc also exotic inflationary models in the 
literature that can achieve tl, > 0, either by violating 
the weak-energy condition [65| or by modifying gravity 
[66l |. Furthermore, although the upper bound = 

applies to inflationary cosmology, it does not apply to 
other cosmological models in which the perturbations 
are produced during a contracting phase {e.g. "pre-big- 
bang" cosmological models, which predict = 
cyclic/ekpyrotic models, which predict ri^ w 2 

It is perhaps worth adding that, instead of consider- 
ing inflation in general terms, one may wish to focus 
on single-field inflation. After all, in 2007, the simplest 
single- field inflation models {e.g. the quadratic inflaton 
potential V{(f>) = {l/2)m^(j)^) continue to agree beau- 
tifully with the current cosmological data sets [sHi, and 
arguably provide the simplest and most compelling avail- 
able explanation of those data sets. Since single-field 
models satisfy the well known "inflationary consistency 



turns out that we can 
' -'^max/S in the upper 



relation" n^{k^^^) -r/8, it 
make the substitution fi^ — 
bound (|2ip . and thereby obtain a somewhat stronger up- 
per bound that is still obeyed by nearly all single-field 
inflationary models. 

To stress that the upper bound on {f) at high fre- 
quencies is very sensitive to the assumed values for w^^y. 
and n^ we plot this upper bound in Fig. [51 for various 
choices of w^^^ and fit^^^^- 

Fig.[2]also shows the bounds and sensitivities from var- 
ious current and future gravitational-wave constraints. 
The LIGO experiment is currently operating at its design 
sensitivity, and has placed an upper bound {f) < 
6.5 X 10~'^ on the stochastic gravitational-wave back- 
ground at frequencies near / ^ 10^ Hz (g^]. The LIGO 
sensitivity is expected to increase by another factor of 
10-100 within the next year or so [63|. Then, within the 
next ten years, Advanced LIGO/VIRGO is expected to 
reach a sensitivity of nl^{f) « 10"^ - 10"^ [13; and 
subsequent generations of ground-based LI experiments 
may do even better. LISA (the first-generation space- 
based LI experiment) is expected to achieve a sensitivity 
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of Hf"{f) « 10~^^ at frequencies near / ^ 10~^ Hz 
[6^ : and BBO (the second-generation space-based LI 
experiment, which is specifically designed to detect a 
stochastic gravitational-wave background) may be able 
to reach a sensitivity of flQ^" {f) w 10~^^ at frequencies 
near / 0.3 Hz [2^, [6^. Pulsar-timing experiments have 
currently placed an upper bound ri§^(/) < 2 x 10~® 
at frequencies between 10~^ and 10~* Hz [l^. In the 
coming years, the Parkes Pulsar Timing Array (PPTA), 
which is already operating, should reach a sensiti vity 
of Qq^^/) w lO"-"^" or better at these frequencies p^ : 
and in the future, the proposed Square Kilometer Array 
(SKA) experiment may improve this sensitivity by an- 
other order of magnitude or more [Hi. Finally, if short- 
wavelength primordial gravitational waves had too much 
energy density, they would spoil the successful predic- 
tions of BBN; so we obtain the standard Bi g B ang Nu- 
cleosynthesis (sBBN) constraint [tI!, Ill III Izl UM, III 
depicted in Fig. [21 and further discussed in Sec. I VIII 



<s= 




f= 100 Hz (LIGO) 

f = 0.3 Hz (BBO) 



- f = 10 3Hz (LISA) 

- f = 10-9 Hz (PULSAR) 



B. Case 2: both r and Q^™ (f) are detected 

If CMB experiments succeed in detecting r, and one 
of the LI (or PT) experiments (at frequency /) also suc- 
ceeds in detecting f2Q*(/), then the master equation ([T|) 
will yield a curve in the {w{f),n^{f)} plane. This is 
illustrated in Fig. [3l In fact, this curve will be slightly 
"fuzzy" due to the non- vanishing error bars on r, CIq" (f), 
Ai, A2, and A3. 

In particular, CMB polarization experiments are ex- 
pected to be sensitive to a tensor-to-scalar ratio as small 
as r = 10~^, or smaller and the sensitivities of 

current and future PT and LI experiments were dis- 
cussed at the end of Sec. IV AI In the top panel of 
Fig. [31 we imagine that ^Iq^^IJ)/!' = 10^^ has been de- 
tected — e.g. r = 0.1 has been detected in the CMB, 
and Oq"(/) — 10^^ has been detected in one of the 
LI/PT experiments — and then we plot the correspond- 
ing constraint curves, assuming that the detection of 
fll^if) occurred at a frequency of lO^ Hz (LIGO), 0.3 Hz 
(BBO), 10-3 Hz (LISA), or IQ-^ Hz (PT). And the bot- 
tom panel of Fig. [3] illustrates the same thing, assum- 
ing that the value of ^lf"{f)/r turns out to be closer 
to the minimum possible value for each experiment: we 
use nl^/r = 10-70.1 = 10-s for LIGO (at 10^ Hz); 
and nl^/r = lO^^VO-l = IQ-^^ for BBO (at 0.3 Hz); 
and nf'/r = 10-11/0.1 = 10-^° for LISA (at 10"^ Hz); 
and nl^/r = 10-"/0.1 = 10-^° for SKA (at IQ-^ Hz). 
Note that, in Fig. [31 the frequency / is different for each 
LI/PT experiment: that is, LIGO places a constraint 
in the {w(/ligo)' "^((/ligo)} pl^'He, while LISA places 
a constraint in the {w(/lisa)i "-iC/LisA)} pla-nc, and so 
forth. 



FIG. 3: Bounds from combining CMB and LI/PT experi- 
ments. (This figure relates to Cases 2, 3, and 4, discussed in 
Sees. IV Bl IV CI and IV Dl l We show examples of the con- 
straints in the {■w{f),nt(f)} plane that follow from CMB 
constraints on r and LI/PT constraints on QQ^{f). In both 
the top and bottom panels, the 4 curves correspond to: 
./ligo = 100 Hz (red dotted); /bbo — 0.3 Hz (green dashed); 
/lisa = 10"^ Hz (grey solid); and /pui^ar = 10"^ Hz (black 
dot-dashed). In the top panel, all 4 curves are plotted as- 
suming QQ"{f)/r = 10"^. So, for example, suppose CMB 
and LI/PT experiments find: (i) f7g™(/) = IQ-** and r = 0.1 
(Case 2); or (ii) Q.1^{f) = 10-* and r < 0.1 (Case 3); or (iii) 
^rif) < 10"** and r = 0.1 (Case 4). Then {w{f),ht{f)} 
must lie: (i) on, (ii) above, or (iii) below the respective curve. 
The bottom panel is similar, but instead of all curves corre- 
sponding to f7Q"(/)/r = 10"^, we take f7Q"(/)/r to be closer 
to the minimum possible value for each experiment: 10"* 
(at /ligo); 10"" (at /^sa); 10"" (at /bbo), and 10"" (at 

/"pulsar ) ■ 

C. Case 3: Qf^{f) is detected, but r is not detected 

In this section, let us suppose that one of the LI/PT ex- 
periments has successfully detected r2o™(/) at some fre- 
quency /; while CMB experiments have only placed an 
upper bound on the tensor-to-scalar ratio: r < r^^^. We 
will mention three possible interpretations of this obser- 
vational situation. 

For the first interpretation, we rewrite the master equa- 
tion dll) as: 

< [^l4^^^4'^^Vmax- (25) 

As in Case 2, this equation defines a curve in the 
{w(f),n^(f)} plane. But, whereas in Case 2 the param- 
eters t«(/) and n^{f) were required to lie on this line, in 
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the present situation the parameters are required to lie 
above the line (see Fig. [3]) . 

For the second interpretation, we rewrite the master 
equation ([T]) as: 



Wif) > ^^mi„(/)' 



where 



and 



^^mi„(/) 



3 V2 



.(/) 



"mi„(/) 



n„(/) = - 



ln[AiA 



nt.„ 
3 



.(/) 



hl[A2 



(26) 



(27) 



(28) 



In other words: if we assume a theoretical upper bound 
for fi-tif), such as the standard inflationary assumption 
^t.max = discussed in Sec. IV A| then we can infer that 
w^f) must exceed the lower bound w,„in(/) given by 
Eqs. dill) and jM]). Furthermore, Eqs. © and ^ allow 
us to infer that the effective equation-of-state parameter 
'w{a) must also satisfy the same lower bound 



wia) > w„in(/) 



(29) 



for some non-empty subset of the range < a < ac- 
And then, if we assume C(a) > (as required by the sec- 
ond law of thermodynamics), we can also infer that the 
ordinary equation-of-state parameter 'w{a) must again 
satisfy the same lower bound 



w{a) > w^iM) 



(30) 



for some non-empty subset of the range ak < a < Gc- 

Eqs. (H?]) and for u),„in(/) as a function of 

^o^(/)/^maxj ^I'S plotted in the top panel of Fig. 21 as- 
suming rnax(/) = (t^c Standard inflationary assump- 
tion, discussed in Sec. IV Ap . The 4 curves correspond 
(from top to bottom) to: FT experiments at / ~ 10~^ Hz 
(black dot-dashed curve); LISA at / 10~^ Hz (grey 
solid curve); BBO at / ~ 0.3 Hz (green dashed curve); 
and LIGO at / ^ 10^ Hz (red dotted curve). Since the 
curves represent w^i^{f)^ the actual value of w{f) must 
lie above these curves. 

For the third interpretation, we rewrite the master 
equation ^ as: 



where 



and 



.i„(/) 



fitif) > rit_„,in(/), 
ln[r^,,^i4-(^V^ir(/)] 



amax(/) = 2 



\3Wmax(/) + 1 



(31) 



(32) 



(33) 



9" 
"O.max 



(f)/r 




3> 
3 



10-15 



10-10 
i^r(f)/rma 



10-5 



FIG. 4: Bounds from combining CMB and LI/PT experi- 
ments. (This figure relates to Cases 3 and 4, discussed in 
Sees. IV CI and IV Dh In both the top and bottom panels, 
the 4 curves correspond to the 4 frequencies: /ligo — 100 Hz 
(red dotted); /bbo = 0-3 (green dashed); /lisa — 10^'* Hz 
(grey solid); and /pulsar ~ 10~^ Hz (black dot-dashed). This 
figure has 2 interpretations. In Case 3, where LI (or PT) 
experiments detect f2g*(/) and CMB experiments obtain an 
upper bound rmax, the "bottom" and "left" axis labels apply, 
and the curves represent i&niin(/) (top panel, with the stan- 
dard inflationary assumption rit.max = 0) and ftt,min(/) (bot- 
tom panel, with the standard primordial-dark-age assumption 
Wmax = 1/3), so the actual values of w{f) and nt[f) lie above 
the curves. In Case 4, where CMB experiments detect r and 
LI (or PT) experiments obtain an upper bound f2o^ax(/)' 
"top" and "right" axis labels apply, and the curves represent 
'Wma.ycif) (top panel, with the standard inflationary assump- 
tion fit ~ 0) and nt,max(/) (bottom panel, with the standard 
primordial-dark-age assumption 'w{f) « 1/3), so the actual 
values of w(/) and nt{f) lie below the curves. 



In other words: if we assume a theoretical upper bound 
for w(/), such as the standard assumption Wj^ax = 1/3 
discussed in Sec. IV Al then we can infer that «((/) must 
exceed the lower bound ^i^^if) given by Eqs. ([32]) and 
([55)) . Furthermore, from Eq. we can infer that the 
actual primordial tensor power spectrum n^{k') must also 
satisfy the same lower bound 



(34) 



for some non-empty subset of the range k^^^^ < k' < k. 

and ([33)) . for fif j^i^if) ^ function of 



Eqs. dsn 

(/)/^inaxj plotted in the bottom panel of Fig. 21 
assuming Wi„ax(/) — 1/3 (a standard assumption about 
the primordial dark age, as discussed in Sec. lV A[) . Again, 
the 4 curves correspond (from top to bottom) to: PT 
experiments at / '--^ 10~^ Hz (black dot-dashed curve); 
LISA at / - 10-3 Hz (grey solid curve); BBO at / ~ 
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0.3 Hz (green dashed curve); and LIGO at / 10^ Hz 
(red dotted curve). Since the curves represent fi^. „!;„(/), 
the actual value of fi^^f) must lie above these curves. 



D. Case 4: r is detected, but flg^{f) is not detected 

Finally, in this section, let us suppose that CMB ex- 
periments have successfully detected a non-zero value for 
r, but LI/PT experiments have only managed to place 
an observational upper bound rig"'(/) < f^o^ax(/) 
frequency /. As in the previous section, we will men- 
tion three possible interpretations of this observational 
situation. 

For the first interpretation, we rewrite the master equa- 
tion III) as 



2 ^3 



(35) 



As in Cases 2 and 3, this equation defines a curve in the 
{w{f),n^{fy} plane. But, whereas in Case 2 the param- 
eters were required to lie on this curve, and in Case 3 the 
parameters were required to lie above this curve, in the 
present case the parameters are required to lie below this 
curve (see Fig. 

For the second interpretation, we rewrite the master 
equation ^ as: 



>(/)< 



where 



and 



Af) 



.(/) 



1/2. 
3 V2 



c(/) 



.(/) 



ln[rAiAg-^^VOg-,,,(/)] 
HA2] 



(36) 



(37) 



(38) 



In other words: if we assume assume a standard value 
for n^{f), then we can infer that w(/) must be less than 
the upper bound w^^^{f) given by Eqs. (|37|) and ((38)) . 
In inflation, the primordial gravitational wave spectrum 
is extremely flat, so that the "standard" value may be 
taken as n^[f) « 0. In fact, the standard inflationary 
gravitational wave spectrum has a slight negative tilt, 
n((fc) = — 2e(fc), but it is small enough that we can ig- 
nore it for the purpose of keeping the present discussion 
simple. It is enough to note that the slight fuzziness in 
the standard inflationary value n^{f) w leads to slight 
fuzziness in the inferred upper bound w^-^^^{f). 

Eqs. ((57|) and (pS)) . for w^^^{f) as a function of 
r2Q^a^x(/)/r, are plotted in the top panel of Fig. 21 as- 
suming the standard inflationary value n^{f) k, 0. The 
4 different curves correspond to the different LI/PT fre- 
quency bands, as already described for Case 3 in Sec. IV Cl 
But, in Case 3, these curves represented w^^^{f), so that 
the actual value of w(/) was required to lie above the 



curves. And now, in Case 4, these same curves represent 
^maxif)^ SO that the actual value of w{f) is required to 
lie below the curves. 

For the third interpretation, we rewrite the master 
equation ([T]) as: 



fltif) < nt.max(/)> 



(39) 



where 



and 



. ln[rAi4^^Vl7§-_(/)] 

"t,max(/) = rm (40) 



K/) = 2 



HA3 

( Mf) - 1 

\3wif) + 1 



(41) 



In other words: if we assume a standard value for w(/), 
then we can infer that n^{f) must be less than the up- 
per bound „ja^(/) given by Eqs. (gO]) and (gl]). The 
most common picture of the post-inflationary universe 
is that, after reheating completes, the universe settles 
quickly into ordinary quasi-adiabatic radiation-like ex- 
pansion [bJ], so the "standard" value may be taken 
as w{f) w 1/3. In fact, even during standard quasi- 
adiabatic radiation-like expansion, various effects — no- 
tably conformal anomalies [gO, [ZO] and the evolution of 
g^, and g^,^ with time [1^ — cause w to drop slightly be- 
low 1/3, but these corrections are usually small enough 
that we can ignore them for the purposes of keeping the 
present discussion simple. It is enough to note that the 
slight fuzziness in the standard value w(/) « 1/3 leads to 
a slight fuzziness in the inferred upper bound fif max(/)- 



Eqs. (gOl) and dH]), for 



i.max 



(/) as a function of 



^o^ax(/)/^' are plotted in the bottom panel of Fig. 21 
assuming a "standard" primordial dark age: w{f) « 1/3. 
The 4 different curves correspond to the different LI/PT 
frequency bands, as before. But, in Case 3, these curves 
represented h^^:^^{f), so that the actual value of n^{f) 
was required to lie above the curves. And now, in Case 4, 
these same curves represent max(/)i that the actual 
value of nt(/) is required to lie below the curves. 



VI. OBSERVATIONAL CONSISTENCY CHECK 

Suppose that a pulsar-timing experiment, or a laser- 
interferometer experiment like LIGO or LISA, detects a 
non-zero value for {f) that is far above the expected 
upper bound ^ 10"^^ which follows from assuming "stan- 
dard" inflation plus a "standard" primordial dark age 
(see Sec. IV Al and Fig. |5]). If we wish to interpret this 
as a detection of the primordial gravitational- wave back- 
ground, then we should expect it to satisfy the following 
rough consistency check. 

If the unexpectedly high value of Vl^ [J) is really due 
to an unexpectedly high value of w{f), or an unexpect- 
edly high value of fitif), or both, then 17q"(/) should 
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be very "blue," i.e. rapidly rising with frequency. This 
point should be intuitively clear from a glance at Fig. [21 
but let us be a bit more quantitative. The standard ex- 
pectation is that all four terms on the right-hand-side 
of Eq. are nearly zero, and hence ilQ^(/) is nearly 
frequency- independent. But if the detected signal is ac- 
tually due to an unexpectedly high value of w{f), then 
the first term dominates the right-hand-side of Eq. (fTU)) . 
and we expect 



bbn 



din/ 



1 



(42) 



where w^^^^ is given by Eqs. p7)) and p8)) . And, similarly, 
if the detected signal is actually due to an unexpectedly 
high value of n^(f), then the second term dominates the 
right-hand-side of Eq. (fTO)) . and we expect: 



(43) 



where ni^min is given by Eqs. ([32]) and ([33]) . These ex- 
pectations can be checked within the frequency band of a 
single experiment, or by comparing two different interfer- 
ometers with two separated frequency bands (like LIGO 
and LISA). 

Note that this is just a consistency check - it does not 
rule out the possibility that the detected gravitational- 
wave signal is produced by some other source, such as a 
cosmological phase transition, cosmic strings, or an unan- 
ticipated astrophysical source. Furthermore, we have 
been careful to use the term "expect" rather than "pre- 
dict" in this section, since it should be clear that Eqs. 

and P5)) are not firm predictions. Nevertheless, they 
arc sufficiently strong expectations that — depending on 
whether or not they are confirmed — they could signif- 
icantly increase or decrease our confidence in the "Case 
2" or "Case 3" interpretations discussed in Sees. IV Bl and 

Ecl 



VII. CMB + BBN CONSTRAINTS 

In this section, let us suppose that CMB experiments 
have succeeded in detecting r, and combine this infor- 
mation with the well-known "standard Bi g B ang Nucle- 
osynthesis" (sBBN) constraint [zl Izl, III, III, IzH, Izll : 



1.8 X lO"^'^ Hz, corresponds to the mode that 
was on the Hubble horizon {ky^^^ = Obbn^bbn) 
time of BBN, while the upper integration limit f^^^^ cor- 
responds to the high-frequency cutoff of the primordial 
gravitational-wave spectrum. For example, if the pri- 
mordial gravitational-wave spectrum was generated by 
inflation, then the spectrum cuts off exponentially fast 
for k > kend, where k^^^ = a^^^H^^^ is the comoving 
wavenumber that was on the Hubble horizon at the end of 
inflation. This corresponds to the present-day frequency 
given by Eq. 

Although; the sBBN constraint (HJ) is technically an 
integral constraint (non-local in frequency space), in 
practice it effectively acts as an algebraic constraint (lo- 
cal in frequency space) of the form fl^"" {f) < 1.5 x 10^^ 
for /jjjjjj < / < /end- ^o^(/) Can only exceed this bound 
by having a very narrow spike with ((S/)//o 1, where 
/o is the peak of the spike, and Sf is its characteristic 
width; but (as far as we are aware) there are no known 
mechanisms for producing such a narrow spike in the pri- 
mordial gravitational-wave spectrum, and we will neglect 
this possibility. 

Thus, for any frequency / in the range f^^^^^ < f < 
/(,„(j, we can directly use all of the equations from "Case 
4" in the previous section, as long as we set f^o^ax(/) = 
1.5 X 10~^ in those equations. Furthermore, to maxi- 
mize the length of the "lever arm" between the CMB 
and BBN constraints, let us consider the case k — > k^^^ 
and / /end- Then Eqs. (j^OI and (j^Tl) define a curve in 
the {w(/e„d),^t(/end)} pl^^c (shown in Fig.©, and the 
actual values of w(/ond) "-((/end) niust lie below this 
curve. 

Note, in particular, that the constraint curve hardly 
varies as r varies over the range of realistic future de- 



tectability 10 



< r < 10" 



Furthcrmore, for r 



{10""'^, 10~^, 10^'^}, respectively: (a) if rij(/) is assumed 
to take its "standard" inflationary value {n^{f) ~ 0, see 
Sees. IV Al and IV C|) . then we obtain the upper bound 
wUend) ~ {0.54,0.57,0.60}; and (b) if w(/) is assumed 
to take its "standard" post-inflationary value {w{f) w 
1/3, again see Sees. IV Al and IV Cp . then we obtain the 
upper bound n^if^^^) < {0.36,0.40,0.43}. 

These results (particularly Fig.[5l) are new, and model- 
independent, constraints on the early universe that will 
take effect as soon as CMB polarization experiments de- 
tect a non-zero value for r. 



^T{f)j < 1.5 X 10-^ 



(44) 



VIII. CONCLUSION 



Note that this constraint only applies to the part of the 
present-day gravitational-wave spectrum that was gen- 
erated prior to Big Bang Nucleosynthesis; and the inte- 
gral only runs over frequencies / corresponding to co- 
moving wavenumbers k that were already "inside the 
Hubble horizon" (fc > Obbn^bbn) time of Big 

Bang Nucleosynthesis (BBN) at photon temperature 
T ~ 1 MeV. In particular, the lower integration limit. 



As far as the early universe is concerned, most peo- 
ple think about upcoming cosmic-microwave-background 
(CMB) "B-niode" polarization experiments with the fol- 
lowing goal in mind: to measure one crucial number, 
r, which physically corresponds to measuring the en- 
ergy density of the universe, roughly 60 e-folds before 
the end of inflation. But these B-mode experiments 
will actually achieve significantly more than this: they 
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FIG. 5: Bound from combining sBBN and CMB constraints. 
If CMB experiments detect r, then the sBBN gravitational- 
wave constraint immediately requires {ti'(/end)i '^tC/end)} to 
lie below the curves shown in the figure. From top to bottom, 
the curves correspond to: r — 10"'^ (black dotted curve), 
r = lO^'^ (purple dot-dashed curve), and r — 10^^ (green 
solid curve). Note that the curves are very insensitive to 
r; they hardly move as r varies over the range in which 
it can be realistically detected by CMB polarization exper- 
iments (10""^ < r < lO"'^). The horizontal and vertical 
dashed lines point out that, for r = {10^^, 10"^, 10~^}, re- 
spectively: (a) if nt{f) is assumed to take its "standard"' 
value (« 0), then w)(/end) ^ {0.54,0.57,0.60}; and (b) if 
w{f) is assumed to take its "standard" value (~ 1/3), then 
nt(/cnd) < {0.36,0.40,0.43}. 



should also be viewed as half of a two-pronged exper- 
iment to detect or constrain the early-universe param- 
eters w{f) and fif^f), as we have described in detail. 
(The other "prong" of this two-pronged experiment is 
a higher-frequency gravitational-wave constraint com- 
ing from laser-interferomctcr experiments, pulsar-timing 
measurements, or standard BBN.) For example, if and 
when CMB experiments detect a non-zero value for r, 
they will immediately obtain a supplementary (and re- 
markably strong) constraint in the {w(/end)' "■t(/end)} 
plane, as shown in Fig. [5l Since quantitative and model- 
independent constraints on the early universe are noto- 
riously hard to obtain, and we only have a handful, the 
possibility of obtaining this "supplementary" constraint 
is exciting. 

We have argued that combining large- wavelength con- 
straints on r (from CMB experiments) with small- 
wavelength bounds on rig^(/) (from LI and PT experi- 
ments, and sBBN constraints) provides the strongest way 
to constrain (or detect) the existence and properties of a 
possible exotic "stiff energy" component (with w > 1/3) 



[13, H m, M, mi, m, S, H H m t^at couU have 
dominated the universe for some period during the pri- 
mordial dark age between the end of inflation and the 
BBN epoch (see Fig.[T]). 

We have derived several useful and general formulae 
for relating primordial gravitational-wave constraints at 
different frequencies, and have shown how these relation- 
ships connect to the uncertain physics of the early uni- 
verse. In Figs, m [3l m and [5l we have shown the con- 
straints that will be placed on the parameters w{f) and 
hi{f) by combining various pairs of gravitational- wave 
constraints, depending on the observational situation 
(that is, depending on whether CMB and/or LI/PT ex- 
periments detect the primordial gravitational- wave back- 
ground, or only place upper limits). 
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APPENDIX A: DERIVING THE MASTER 
EQUATION 

The goal of this Appendix is to derive Eq. (|A34p . which 
is equivalent to Eq. ([1]), the master equation upon which 
most of the analysis in this paper relies. As a useful inter- 
mediate result, we also obtain Eq. (|A33p . an expression 
for the tensor transfer function r^ (fc). 

The derivation is broken into 4 parts. In the first part, 
we review some background material about cosmologi- 
cal gravitational waves, leading to the presentation of 
Eq. (|A7[) . The second and third parts are devoted to 
rewriting the factors I)j(fc) and A^(fc,ri), which appear 
in Eq. (|A7[) . Finally, in the fourth part, wc collect and 
summarize our results in Eqs. (|M3)) and jMH- 

1. Background material 

Let us start by introducing some notation, and re- 
viewing some basic facts about cosmological gravitational 
waves (tensor perturbations [ttI). For more details, see 
Sec. 2 in Ref. [G^. 

Tensor metric perturbations in a spatially-flat FLRW 
universe are described by the line element: 

ds^ = a2(r) [-dr^ + (% + % (x, r))dx'dx^] , (Al) 

where x is a comoving spatial coordinate, r is a conformal 
time coordinate, a{T) is the FLRW scale factor, and the 
metric perturbation hij{x., t) is transverse (hijj = 0) and 
traceless (ha =0). In this Appendix we follow the con- 
vention that repeated indices (i or j) are summed from 
1 to 3. 
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The tensor power spectrum A^(fc, r), which represents 
the contribution by modes of comoving wavenumber k 
to the expectation value {hij{x.,T)hij{x.^T))^ is defined 
through the equation 



(^■(x,T)/ly (x,t)) 



dk 



(A2) 



Note that the expectation value of the left-hand-side is 
actually independent of x, since a perturbed FLRW uni- 
verse is statistically homogeneous. 

CMB and LI experiments measure A^(fc,r) at very 
different comoving wavenumbers k and very different 
conformal times r. In particular, whereas LI experi- 
ments measure the present-day tensor power spectrum 
A^(fc, To), CMB experiments may be thought of as mea- 
suring the primordial tensor power spectrum A|(fc,ri). 
(Here tq denotes the present time, and r,; denotes a very 
early time, before any modes k of interest have had a 
chance to re-enter the Hubble horizon.) And whereas LI 
experiments are sensitive to high comoving wavenumbers 
(corresponding to length scales smaller than the Solar 
System), CMB experiments are sensitive to low comov- 
ing wavenumbers (corresponding to large length scales, 
comparable to the present-day Hubble radius). CMB 
constraints on the primordial scalar and tensor power 
spectra, A^(fc,Ti) and A|(fc,T,;), are usually quoted at a 
fiducial "pivot" wavenumber /ccmb m the CMB waveband. 
For example, the WMAP experiment uses k^^^-^^/ao = 
0.002 Mpc~^, where ao is the present-day (t = tq) value 
of the FLRW scale factor. 

Although it is often convenient, from a theoretical 
perspective, to work with the tensor power spectrum 
A^(A:,r), LI experiments usually report their results in 
terms of the present-day (r = Tq) gravitational- wave en- 
ergy spectrum: 



where 



f 



1 dp^ 



1 k 
2n oq 



(A3) 



(A4) 



is the present-day physical frequency of a gravitational 
wave corresponding to the comoving wavenumber k. 
Note that the present-day energy spectrum il,^'"{f) is 
related to the present-day power spectrum Af^{k,To) 
through the equation 
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27r/ 



A',(fc,ro) 



(A5) 



(see Sec. 2 in Ref. [g^l for a detailed derivation). 

The present-day tensor power spectrum A^(fc,ro) is 
related to the primordial tensor power spectrum A^(fc, t^) 
via the relation 



A2(fc,ro)=r„(fc)A2(fc,r,), 



(A6) 



where this equation defines the "tensor transfer function" 
r^(fc). Combining Eqs. (|X5)) and (K6\i we obtain 



'if) = — 



27r/ 



T^{k)Al{k,n). 



(A7) 



This is the master equation describing the present-day 
gravitational- wave energy spectrum ilQ^(/) on LI scales. 
The rest of this section is devoted to rewriting this equa- 
tion in a more concrete and useful form. In the next 
two sections, we re-express the two factors T)i(fc) and 
Af^{k,Ti), respectively. 



2. Rewriting the factor T,-^{k) 

First, let us focus on rewriting the factor Tf^(k). In 
this paper, we make use of the general expression for the 
tensor transfer function Ty^{k) derived in Ref. [B^I- As 
explained in Ref. [i^l , the tensor transfer function (fc) 
may be factored into the form 



r,(fc) = ^Ci(fc)C2(fc)C3(fc). 



(A8) 



The overall factor of 1/2 comes from averaging over the 
oscillatory factor cos^(fcr + 4>{k)) which appears in the 
tensor transfer function but is unresolvable in any forsee- 
able LI experiment [zl] . Each of the remaining 3 factors 
{Ci(fc), C2(fc), C3(fc)} has a simple physical meaning and 
is derived in detail in Ref. (60j . Here we just quote a few 
key results. 

As we shall see, the expression (|A8[) for Tj^{k) is dom- 
inated by the factor Ci(fc) ^ 1, while the other two 
factors, 6*2 (A:) and C3(A:), represent modest 0(1) correc- 
tions. 

The factor Ci (fc) is given by 



Ci(fc) = 



1 



(A9) 



where z^, is the redshift at which the mode k re-enters 
the Hubble horizon (fc = aH) after inflation. We shall 
return to this factor below. 
The factor C2{k) is given by: 



C2ik) 



r'(«fe + 1/2) 



2 

a I. 



(AlO) 



where T{x) is the gamma function, and we have defined 

2 



a,. 



1 + 3w,. 



(All) 



Here w^. is the effective equation-of-state parameter at 
redshift z^,, and is given by 



Wi. = Wi. 



STrG^Cfe 



(A12) 
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FIG. 6: The correction factor C2 as a function of the effective FIG. 7: The correction factor C3 as a function of the free- 
equation-of-state parameter w, as given by Eqs. (|A10|I and streaming fraction = pI" I Pk'^ , as given by Eqs. (|A13|) 
(|A11|) . Note that C2 = 1 when m) = 1/3. and (fAT4|l . 



where w^. = p^j/pj, is the usual equation-of-state param- 



eter (i.e. the ratio of the total cosmological pressure 
to the total cosmological energy density p^.), H^. is the 
Hubble expansion rate, Cfc is the bulk viscosity of the 
cosmological fluid (see Sees. 2.11 and 15.11 in Ref. [6l| ) 
— and; as their subscripts indicate, all of these quantities 
are evaluated at redshift z^. C2 is plotted in Fig. [HI Note 
that the expression (jA10[) for C2{k) is valid as long as 
the effective equation-of-state parameter w is not chang- 
ing rapidly relative to the instantaneous Hubble time at 
redshift z^.; see Ref. [gO] for the meaning of C2(fc) more 
generally. 

The factor C^^k) captures the modification of the 
primordial gravitational-wave signal due to tensor 
anisotropic stress Tr^j (e.5. from free-streaming relativis- 
tic particles) in the early universe. In particular, in the 
important case that the effective equation-of-state near 
Zk is radiation- like [w « 1/3), free-streaming relativistic 
particles with energy density p^'^ damp the gravitational 
wave spectrum by the factor: 



C3(fc) = A\k) 



where 

A{ky. 



10 (98f73^-589rj2 



9380^2^^-55125) 



7 (15-l-4f}^J(50 + 417^.J(105-h417^J 



and fijj, = Pfe^/p™* is the fraction of the critical density 
that is relativistically free-streaming at redshift z^,. C3 is 
plotted in Fig. [71 

In the remainder of this section, we focus on obtaining 
a more explicit expression for (1 -I- z^.), and hence for 



(A13) 



(A14) 



Ci{k). To do this, let us proceed carefully as follows: 



"0^0 



(1 + Po 



crit 



mat rad „crit 
PeqPc Pk 



/I I ~ ^ 2 ^crit ^rad ^crit ' 
\^ ^ '^kl Po Peq Pc 



(Al5a) 
(A15b) 
(Al5c) 



where p" 



„rad 



and p™^ denote the matter density. 



radiation density, and critical density, respectively. In 
the first line (jA15ap . wc have used the fact that k — 
^k^k by definition. In the second line (jA15bp . we have 
used the definition of the redshift z to write a^/ai^ = 
(1 -I- z^,), and the definition of the critical density p'^"* to 



write Hl/Hl = p-'Vpg- 



9o"*. In the third line (|A15c|) . we 
have used the fact that p™''* = p^^^ at the time Teq of 
matter-radiation equality; plus the fact that the universe 



„rad 



is radiation dominated at so that p"' 

We have introduced the time Tc to parametrize our 
threshold of ignorance: it represents the earliest time at 
which we know that the universe was already radiation 
dominated. But, for all we know, the universe prior to 
Tc may not have been radiation-dominated: e.g. an ex- 
otic "stiff" component with w > 1/3 may have domi- 
nated the cosmic energy budget. The present agreement 
between the theoretical and observational understand- 
ing of Big Bang Nucleosynthesis (BBN) strongly suggests 
that the universe was already radiation dominated dur- 
ing BBN (i.e. at the time t^\^^ when the temperature 
was T = 1 MeV), so it currently makes sense to choose 
Tc ~ ''bbn- But, for the sake of generality, we leave Tc un- 
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fixed in this section, since we can easily imagine future 
developments — such as an improved understanding of 
primordial baryogenesis — that would make an earlier 
time Tc <C "T^hn ^ more appropriate choice. Note that for 
the wavenumbers of interest in this paper — e.g. those 
measured by laser interferometer experiments — the tem- 
poral ordering is Tk < Tc < Teq < Tq. 

Three density ratios appear on the right-hand side of 
Eq. (|A15c|) . Let us rewrite each of these in turn. The 
first density ratio is trivially rewritten as: 



Peg 



(A16) 



where fi™''* 



„mat / ^crit 

Po / Po ■ 



In order to rewrite the second density ratio appearing 
in (jAlScp . let us pause for a moment to review a few 
properties of an expanding bath of radiation. The ra- 
diation bath has energy density p and entropy density s 
given by (see Sees. 3.3 and 3.4 in Ref. [6^, and especially 
Eqs. (3.61), (3.72)): 

Piz) = ^7t'9Az)THz), (Al7a) 
K^) = ^Tr'g.MTHz), (A17b) 



To rewrite the third density ratio appearing in 
Eq. (|A15c[) . note that conservation of stress-energy 
(T'^".^ = 0) in the early universe {i.e. in a spatially-flat 
FLRW universe) implies the continuity equation: 



dp"" 



p- 



, ..da 
^3[l-f u.(a)] — . 



(A20) 



Here w(ci) is the effective equation-of-state parameter: 

87rGjvC(a) 



w{a) = w{a) 



H{a) 



(A21) 



where w{a) = p{a)/p{a) is the usual equation-of-state 
parameter [i.e. the ratio of the total cosmological pres- 
sure p{a) to the total cosmological energy density p(a) ~ 
p'^"^{a)], H{a) is the Hubble expansion rate, and C,{a) 
is the bulk viscosity of the cosmological fluid (see Sees. 
2.11 and 15.11 in Ref. [6l|). Integrating Eq. ([X20l) from 
flc = «('rc) to flj^ = a(T^.) yields 



Pjr 

^crit 



exp ■ 



3[1 



(A22) 



Alternatively, we can deflne an averaged effective 
equation-of-state parameter w(/) through the equation 



where T(z) is the temperature at redshift z. These equa- 
tions may be taken as the definition of the quantities 
g^{z) and which represent the effective number of 

relativistic degrees of freedom in the radiation bath at 
redshift z, as measured by the energy density p{z) or the 
entropy density s{z), respectively. If the radiation ex- 
pands quasi-adiabatically — as is usually the case in the 
early universe — then the entropy a'^s remains constant 
(to a very good approximation). When this is true, Eqs. 
(|A17ap and (jA17b|) imply that the energy density of the 
radiation bath redshifts as 



p^^'^jz,) ^ g,{z,) gt's\z^) 
9M)gti\z^) 



l + Zl 
1 + 22 



(A18) 



In particular, since the standard radiation epoch begins 
prior to z^, the radiation expanded quasi-adiabatically 
during the epoch > z > z^q, so we can write: 



„rad 



f \ 4/3/ N 



eq 



9*{Zeq) gti^{z^ ) 



1+ 
1 + 2, 



eq 



(A19) 



Note that if we know the phase-space distribution func- 
tions describing each particle species in the radiation 
bath, then we can compute the quantities g^{z) and 
9*s{z) directly — again see Sees. 3.3 and 3.4 in Ref. [g^I 
for more details. For example, if all relevant particle 
species are in thermal with one another at temperature 
T, then g^ = + {7/8)Nf and 5„ ^ + {7/8)Nf, 
where 7V{, and Nf are the total number of relativistic 
bosonic and fermionic degrees of freedom, respectively. 



P%"' 



1 + 

1 + 2, 



3[l+u.(/)] 



(A23) 



Comparing Eqs. (|A22|) and (|A23|) . we see that w{f) is 
the logarithmic average of the effective equation-of-state 
parameter w{a) over the range a^, < a < a^: 

^^^) = rTV^ r^D(a)-. (A24) 

Note that if wia) is an a-independent constant over this 
range (which we DO NOT assume in this paper) then it 
becomes equal to w. 

Finally, we can plug Eqs. (|A16p . (|A19j) . (|A23P into the 
right-hand-side of Eq. (|Al5cp . solve for (1 + 2^,), and thus 
find 



cm 



1 



(1 



-"^ 27r/ 



(1 + 2j H, 



-4/(l+3to(/)) 



where we have defined 



(")mat 

"o 



9*{Zc) gti^iZeq) 



^+Zeq 9AZeq) gt'J'^Zc) 



(A25) 



(A26) 



3. Rewriting the factor A^(fc,ri) 

Now let us focus on rewriting the factor ^\{k, Ti), i.e. 
the primordial tensor power spectrum on short wave- 
lengths. 
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Recall that the tensor spectral index n^{k) is defined 
as the logarithmic slope of the primordial tensor power 
spectrum A^(fc, r^): 



Finally, it is conventional (and also convenient, for 
certain purposes) to define the tensor-to-scalar ratio r 
through the equation: 



d[lnk] 

Integrating this equation, we obtain 



dk' 



(A27) 



(A28) 



cmb' i 



(A31) 



Combining Eqs. (|M9)) and we can rewrite 

A^ (fc, Ti), the primordial tensor power spectrum on short 
wavelengths, in the form: 



where Al^{k^^y^,Ti) is the primordial tensor power spec- 
trum, evaluated at the CMB "pivot" wavenumber k^^^^. 
Alternatively, we can define an averaged spectral index 



through the equation 



(A29) 



Al{k,n) = r A^(fc„„b,T.)[fc/fcc„J"'^^^- (A32) 



4. Recapitulation 



In other words, the effective spectral index n^{f) is noth- 
ing but the logarithmic average of the actual spectral 
index nf.{k) over the wavenumber range from k^^^^ to fc: 



Now let us assemble our results. Plugging Eq. (jA25|) 
into the right-hand-side of Eq. (|A8p . wc obtain the tensor 
transfer function T/^ik) in the useful form: 



ntif) 



ln[fc/fc, 



cmbJ 



dk' 



(A30) 



Note that if is a fc-independent constant over this range 
(which we DO NOT assume in this paper) then it be- 
comes equal to fi^. 



C2ik)C,ik) 
2(l + zJ2 



(l + z,) H, 



-4/(l-H3u.(/)) 



(A33) 



Then, plugging Eqs. (|A32p and (|A33P into Eq. ([All), we 
obtain our final result: 



J 



nS^(f) ^ ^A|(fccmb.-^»)g2(fc)C'3(fc)7 

U) 24 



(1 + Zc) Ho 



"(/) 



fhif) 



(A34) 



r 



which is equivalent to Eq. ([T]) in the text. 



In Eqs_^33| a nd (|X34l) : C2{k) is given by Eqs. 
(lUOlj ([XTT|l . and (Uni; Csik) is given by Eqs. 
and (jA14[) : w{f) is given by Eq. (IA24P : 7 is given by 
Eq. (|A26|) : n^{f) is given by Eq. (jA30|) : and we have 
defined 



(A35) 



Note that, if the quantities C2{k), Gz{k), w{f) and 
n^{f) are only weakly fc-dependent, then the frequency- 
dependences of the tensor transfer function Tf^[k) and the 
present-day gravitational- wave energy spectrum VL^ [f) 
are given roughly by T,^lk) oc (27r//Ho)-''/(i+3'^(^» and 
nl""{f) oc {2nf/Ho)^^f'^+^tU)^ respectively. 



APPENDIX B: DERIVING THE FREQUENCIES 

/c AND /end 

In this appendix we derive Eqs. ([T^ and for 
and /end, respectively. Let us start with Eq. ([12]) for fc- 
We start by writing 



(27r/.)2 



,crit 



(1 + Zc? 



Hi 



mat rad 



(1 + zc)2 p-' 



rad 



(Bla) 
(Bib) 
(Blc) 



If these steps are unclear, see Eqs. (|A15ap . (jAl5b|) . 
(|A15c[) and the paragraph that follows them. Now, sub- 
stituting Eqs. (jAiep and (|A19p into Eq. (jBlc^ . and solv- 
ing for fc, we obtain Eq. ([T^ as desired. 

Next let us derive Eq. ([T5)) for f^^^- To begin, note 
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that we can write Af^{k^^^,Ti) in two different ways. On 
the one hand, using Eqs. (|A29p and (|A3ip . we can write 



end 



cmb 



(B2a) 



On the other hand, we can use the well-known inflation- 
ary formula 



A^(fc,„d,T.) = 647rG 



N 



H. 



end 



2tt 



(B2b) 



where our conventions match those of the WMAP ex- 
periment {e.g. see Eq. (A13) in Ref. [l^l)- Comparing 
Eqs. (|B2a|) and ((B2bl) . we see that 



rr2 _ fl^MKn^ f ^on6 
-"end — 



IGttG 



N 



cmb 



fit(/o„d) 



Next, from the definition of /„„hi wc can write 



end ' 
2 



a 



end/ 



end j:r2 
,2 end- 



a, 



(B3) 



(B4) 







Note that, since the first factor on the right-hand-side of 
Eq. ((B4)l is nothing but Ci{k^^^^), we can use Eq. (|A25|) 
to rewrite it as: 



^end 



1 



end 



(l + zj 



-4/(l+3tf.(/„„J) 



(B5) 



Finally, we can plug Eqs. (|B3p and (jBSp into the right- 
hand-side of Eq. (|B4p . and solve for /j,^^^ to obtain 
Eq. (fT3|) as desired. 



APPENDIX C: NUMERICAL FORMULAE 

This appendix lists a few results that are useful for 
converting our algebraic expressions into numerical re- 
sults and plots. 

At matter-radiation equality, we have the standard 
values g^{zeq) = 3.3626 and g^s{zeq) = 3.9091; and 
at BBN (when the temperature is T = 1 MeV) we 
have the standard values 5*(zbbn) = 5*s('^bbn) =^ 10.75 
[63l |. From WMAP 3rd year data alone, we know 
/^niKuJ = (2-04 ± 0.14) X 10"^ and from WMAP 
3rd year data plus lyman-alpha-forest data we know 
^^(^cmb) = (2-24 ± 0.11) X 10-9 [Zl, where these val- 
ues for A.^(A:^j^,^) are quoted at k^^^^ ~ 0.05 Mpc"^. 
The present day value of the Hubble expansion rate 
may be written as ifp = (3.24)/i x 10~^®Hz, where the 
Hubble parameter h « 0.72 is a fudge factor that ab- 
sorbs the uncertainty in the measurement of Hq. Thus 
we can write (27r//i?o) = (1.94//i) x W^{f /YLz). If 
/c^^b = 0-002 Mpc-\ then (fc^^^/ao^^o) = 6.00/ft.; and 
if ^emb = 0.05 Mpc-\ then {k^^JaoHo) = 150.0//i. 
The redshift of matter-radiation equality may be writ- 
ten as (1 + z^g) = 2.3 X lO-^l^lf^'/i^; and the redshift 
of BBN (when the temperature was T = 1 MeV) is 
(l + 2,,J = 5.9xl09. 
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